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We consider the zero-temperature random-field Ising model in the presence of an external field, on ladders and 
in one dimension with finite range interactions, for unbounded continuous distributions of random fields, and 
show that there is no jump discontinuity in the magnetizations for any quasi-one dimensional model. We show 
that the evolution of the system at an external field can be described by a stochastic matrix and the magnetization 
can be obtained using the eigenvector of the matrix corresponding to the eigenvalue one, which is continuous 
and differentiable function of the external field. 
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The nonequilibrium random field Ising model (RFIM) was 
proposed by Sethna et al. 1 as a model for hysteresis and 
Barkhausen noise in ferromagnets. Since then, there has been 
a considerable theoretical interest in the nonequilibrium re- 
sponse of the model. I 2 - 3 - 4 - 5 ' 6 ' 7 ' 8 Sethna et al. showed that in 
the mean-field limit, if the strength a of the quenched random 
field is large, the average magnetization per site is a contin- 
uous function of the external field, but for small a, it shows 
a discontinuous jump as the external field is increased. How- 
ever, in the mean-field limit, there is no hysteresis above the 
critical disorder a c , i.e., the magnetization follows the same 
curve in the increasing and decreasing field. This shortcom- 
ing of the mean-field limit can be overcome by working on 
a Bethe lattice, where the the model is solved exactly for the 
hysteresis 2 and avalanche size distribution. 3 Interestingly, the 
behavior of hysteresis loops on a Bethe lattice depends non- 
trivially on the coordination number z, so long as the dis- 
tributions of the random fields are unbounded continuous. 
For z — 3 the hysteresis loops show no jump discontinuity 
of magnetization even in the limit of small disorder, but for 
higher z they do. This z dependence of the hysteresis persists 
even for euclidean lattices. 4 The natural question to ask is to 
whether the existence jump discontinuity in magnetizations at 
low disorder depend only on the coordination number or also 
on the dimensionality of the lattice. In this paper we answer 
this question by showing the nonexistence of jump disconti- 
nuity in magnetizations for any quasi-one-dimensional RFIM, 
irrespective of the coordination number, for unbounded con- 
tinuous distributions of random fields. 

We first consider a two-leg ladder (shown in Fig. Q of 
length N. At each vertex there is a Ising spin s* = ±1 which 
interacts with nearest neighbors through a ferromagnetic in- 
teraction J, and coupled to the on-site quenched random field 
hi and the homogeneous external field h. The random fields 
{hi} are drawn independently from a unbounded continuous 
distribution (f)(hi). The Hamiltonian of the system is 



H 



-J 2J SjSj - ^ hjSj - h Si 



(1) 



<%J> 



The system evolves under the zero-temperature Glauber 
single-spin-flip dynamics: 9 a spin flip is allowed only if it low- 
ers the energy. We assume that the rate of spin flips is much 



larger than the rate at which h is changed, so that all flippable 
spins may be said to relax instantly, and any spin s, always 
remains parallel to the net local field li at the vertex 



Si = sign(£i) = sign(J^ Sj + hi + h). 



(2) 



We start with h = — oo, when all the spins are down, and 
slowly increase the field. We are interested in the average 
magnetization of the system as a function of external field, 
away from the two ends of the ladder in the thermodynamic 
limit of N — > oo. 

Under the single spin-flip dynamics, for ferromagnetic cou- 
pling (J > 0), the system exhibits return point memory: 
if we start with with the state where all spins are down (at 
h = — oo), then state of the system at an increased field h(T) 
at a later time T, does not depends on the detailed time de- 
pendent of the external field hit), and is same for all histories 
as long as the condition h(t) < h(T) for all earlier times is 
obeyed. 1 Therefore, to find the magnetization at h, we start 
with h = — oo, so that all spins in the system are down and 
then increase the external field to h in a single step. At that 
field several spins would become unstable. But for J > 0, if 
we start with any stable configuration, and then increase the 
external field and allow the system to relax, the final stable 
configuration reached is independent of the order in which the 
unstable spins are flipped, as flipping of a spin can only in- 
crease the local field at its neighboring vertices and also in the 
relaxation process no spin flips more than once. Because of 
this abelian property of the spin flip, when the external field 
is increased from — oo to h in a single step, we may choose to 
relax the spins from the two ends of the ladder as follows: We 
first relax the spins at the boundary ends (level 1). Then we 
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FIG. 1: Two-leg ladder 
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relax spins next to the boundaries (at level 2) and so on with 
increasing levels as we move towards the center of the lad- 
der. In the relaxation process, if a spin at level x flips up, we 
check all the spins at lower levels for possible upward flips, 
before flipping the second spin at level x. Note that the state 
of the spins at a level r from the right end is independent of 
the state at level I from the left end, so long as all the spins 
in between levels are kept down. Therefore we can relax the 
spins from the two boundary ends on the left and right halves 
of the ladder independently (which are identical). 

We choose our coordinates on each halves of the ladder 
such that the level x consists of vertices (x,l) and (x,2). 
Let P£(si, S2) be the conditional probability that, in the ear- 
lier spin relaxation scheme, the spins at level x reach a fi- 
nal state {si, S2}, i.e., the spin at the vertex (x, 1) — > si and 
the spin at the vertex (x,2) — > S2', given that the spins at 
level x + 1 are kept down and all the spins at the lower lev- 
els are relaxed and the external field is h. Corresponding to 
four possible final states we get four probabilities Pfi (1, 1), 
P£(1,-1),P£(-1, 1) andP£(-l, -1), which depend on the 
states at level x — 1 and the random fields at level x. Note that 
P^(l, — 1) = P% (— 1) 1)> due to the symmetry of the lattice. 

To illustrate how to obtain the recursions for the proba- 
bilities P?(si,S2)'s, we consider an example where the fi- 
nal {1,1} state at level x is achieved from its initial state 
{ — 1,-1}, when state at level a;— 1 is {1,1}. We denote p m (h) 
with < m < 3 as the conditional probability that the local 
field at any vertex i will be large enough so that it will flip up, 
if m of its neighbors are up, when the uniform external field 
is h. Clearly, for a given distribution of random fields <p(hi): 



p m {h) = / 4>(hi) dh h < m < 3. (3) 

J (3-2ro)J-fc 

Since the spin at vertex (x, 1) has one up neighbor, it will flip 
up with probability pi(h). Now the spin at vertex (x, 2) has 
two up neighbors, so it flips up with probability P2(h). But 
the local field at vertex (x, 1) may not be large enough for the 
spin to flip up when it has only one up neighbor. In this case, 
the spin at vertex (x, 2) flips up first, with probability P\(K), 
and then the spin at vertex (x, 1) flips up, when is has two 
up neighbors, with probability [p2{h) — p\(h)]. Therefore, 
the total probability of the spins at level x flipping up via this 
process is Pf _1) (l, l){pi(h)p 2 (h) + [p 2 (h) - Pl (h)} Pl {h)}. 

The case where the state at level x — 1 is { — 1, — 1}, the flip- 
ping of the spin at vertex (x, 1) might causes the spin at vertex 
(x — 1, 1) to flip up and as a result the spin at vertex (x — 1, 2) 
might flip up, while the spin at vertex (x, 2) is still kept down. 

(x—X) 

We denote the probability of this spin flip process by Q h ' 
If we consider [P% (1, 1), 2P£(-1, 1), P£(-l, -1)-Q%, Q%] 
as a column vector V x (h), then the recursion relations for 
these probabilities can be represented in matrix form 



V x {h) = W h V x -i{h), 



(4) 



where Wh is a 4 x 4 matrix whose elements Wy's are polyno- 
mials in p m (h)'s. 

For a given probability distribution of random fields <fi(hi) 
and value of external field h, we determine p m (h), and then, 



using Eq. we can recursively determine the probabilities 
represented by the column vector V x (K). Note that the matrix 
Wh is column stochasticpS i.e., Wij > and J2i=i w ij = !■ 
Therefore, for large x, the vector V x (h) tends to a limiting 
vector V*(h), which is the eigenvector of the matrix Wh, cor- 
responding to the eigenvalue one (the maximal eigenvalue). 
The variation of V*(h) with respect to the external field h is 11 



dP*(h) =A #dWh 
dh h dh 



V*(h), (5) 



where is the group inverse of the matrix Ah , = I — Wh, 
letting I denote identity matrix. The entries of Af are con- 
tinuous functions of h for continuous distributions of random 
fields. Since according to Eq.|5] the elements of Wh are con- 
tinuous and differentiable (with continuous first derivative), 
V*{h) is a continuous and differentiable (with continuous first 
derivative) function of h. 

To calculate the magnetization at the center vertex (0,1), in 
the limit N — > 00, we keep the spins at level O down and relax 
the full system. The states of the spins at the adjacent levels on 
both sides are given independently by the probability vector 
V*{h). Now we relax the spins at at level O. Let P°(-l, -1) 
be the probability that in this relaxation process the spins at 
level O do not flip up and PP(— 1, 1) be the probability that 
only the spin at vertex (O, 2) flips up. The magnetization at 
the vertex (0, 1) is 



M (h) = 1-2 [P° (-1, -1) + P° (-1, 1)] 



(6) 



In the limit N — > 00, all the vertices deep inside the ladder 
are equivalent. Therefore, Mo(h) gives the average magne- 
tization on the ladder far away from the boundary. Since for 
continuous distributions of random fields V*(h) is continuous 
and differentiable, the average magnetization is also continu- 
ous and differentiable function of h. 

We have compared the theoretical calculation with numer- 
ical simulation results. The points in Fig. [2] show the results 
of simulations for a single realization of quenched Gaussian 
random fields with mean and standard deviation a = 0.5 J 
and a = 1.5J, respectively. We use periodic boundary con- 
dition along the length of a ladder with TV = 2 23 . Different 
runs using different realizations of quenched field give results 
which are indistinguishable at the scale of the graph. The solid 
curves the Fig.|5]are obtained using theoretical calculation. 

The analysis of a two-leg ladder can be extended to a case of 
an 71-leg ladder (with finite n), where the relaxation scheme is 
such that, all the n spins at the same level x are relaxed before 
relaxing the spins at level x + 1. The relaxation process can 
be represented by the probabilities Pfi(s\, S2, . . . , s n ) of the 
states of spins at level x and the probabilities Qf(h)'s for in- 
termediate relaxations. Now these probabilities at level x can 
be obtained from the probabilities at level x — 1, and using 
proper linear combinations these can be expressed by a ma- 
trix recursion relation like Eq. 0. It is simpler to consider the 
n-leg ladder to wrapped around a cylinder (of course the re- 
sultant magnetization is quantitatively different from the case 
of an n-leg ladder, but the qualitative behavior remains un- 
changed), since then many of the probabilities become equal, 
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FIG. 2: Magnetization curves for a two-ladder in the increasing field. 
N = 2 23 . 



due to the rotational symmetry of the cylinder. As an exam- 
ple, consider the case n — 3. The state at a level x, when 
the spins at level x + 1 are kept down, can be represented 
by only using P£(l, 1, 1), P£(l, 1, -1), P£(l, -1, -1), and 
P/f(— 1, —1, —1). The Qi(hy& needed to describe the relax- 
ation at level x — 1, when relaxing the spins at level x are 
obtained as follows: 

(a) Suppose the state at level x — 1 is {1,-1,-1} and the 
spin at vertex (x, 2) flips up and the spin at vertex (x, 3) is 
kept down. We define Qj ^ (/i) to be the probability that the 
spin at vertex (x — 1, 2) flips up and subsequently the spin at 
vertex (x — 1, 3) flips up as a result, (b) Suppose the state at 
level x — 1 is {— 1, —1, —1} and the spin at vertex (x, 1) flips 
up while the spins at vertices (x, 2) and (x, 3) are kept down. 
Now we define: (i) (h), the probability that the spin at 

vertex (x — 1,1) flips up and as a result the spins at vertices 

(x— 1, 2) and (x — 1, 3) flip up. (ii)Q^ (h), the probability 
that the spin at vertex (x — 1 , 1 ) flips up and as a result the spin 
at vertices (x — 1, 2) flips up but the spin at vertex {x — 1,3) 
remains down, (iii) ^(h), the probability that the spin at 
vertex (x — 1, 1) flips up but the spins at vertices (x — 1, 2) 

and (x — 1,3) remain down, and (iv) , the probability 

that the spins at vertices {x — 1, 2) and (x — 1, 3) flip up after 
we flip the spin at vertex (x, 2) and the spin at vertex (x, 3) is 
still kept down. 

Now the elements of the probability vector V x {h) 
are P£(l,l,l), 3P£(1,1,-1), 3[P£(1, -1, -1) - Qf(h)] t 
3Qf (A), P,f (-1, -1, -l)-Qf(/i)-Qf(/i)-Qf(/0, QiW, 
Q^(h), Q%(h) —Qf{h) and Q^(h). The recursion relation for 
V x {h) is represented by Eq.|4]with a 9 x 9 stochastic matrix 
W/j whose elements can be expressed in terms of p m (/i)'s. 
Finally, the magnetization can be obtained using the limiting 
vector V*(h) andp m (/i)'s. 




(i + l.l) 



FIG. 3: (a) Linear chain showing nearest and next neighbor interac- 
tions and (b) graph on which nearest neighbor interaction is equiva- 
lent to the linear chain with nearest and next neighbor interactions. 



We next consider the RFIM in one dimension (ID), with 
"finite range interactions" — each spin interacts with all the 
spins up to a distance n, through ferromagnetic interaction 
J. We take the length of the chain N (later we will take the 
limit N — ► oo) to be a multiple of n. Now we can group n 
spins together in a single level and relax the spins from the 
two boundaries such that we relax all the spins at the lower 
levels before relaxing the spins at an higher level. 

As a concrete example, we calculate the magnetization for 
n = 2, i.e., each spin interacts with the nearest and next near- 
est neighbors. The interactions are shown in Fig. EJa). This 
model is also equivalent to the RFIM with nearest neighbor 
interaction on the graph shown in Fig. |3Jb). Let P^(si, S2) 
be the conditional probability that the spins at level x reach 
a final state {si, S2}, i.e., the spin at the vertex (x, 1) — » s% 
and the spin at the vertex (x, 2) — » S2; given that the spins at 
level x + 1 are kept down and all the spins at the lower levels 
are relaxed. We need to define four more probabilities for the 
further relaxation at level x — 1, as a result of a spin flip at 
level x. 

Consider the following cases of relaxing the spins at level x 
from its original state { — 1,-1}: (a) Suppose the state at level 
x — 1 is { — 1, 1} and the spin at vertex (x, 2) flips up, given 
that spins at vertex (x, 1) and level x+ 1 are kept down. It may 
cause the spin at vertex {x — 1, 1) to flip up and we denote the 
probability of it to flip up by Q*~ l (h). The spin at vertex (x — 
1, 1) remains down with probability P^ _1 (-l, l)-Qf _1 (/i). 
(b) Suppose the state at level x— 1 is {—1,-1} and the spin at 
vertex (x, 2) flips up, given that spins at vertex (x, 1) and level 
x + 1 are kept down. This may finally cause the spin at vertex 
[x — 1, 1) to flip up. The spin at vertex (x — 1, 2) may or may 
not flip up during the relaxation. We denote the probability 
of the spin at vertex [x — 1, 1) flipping up by Q2~ l (h). It 
remains down with probability P^ _1 (-l, -1)— Q^^Qi). (c) 
Suppose the state at level x — 1 is {—1,-1} and the spin at 
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FIG. 4: Magnetization curve in the increasing field for linear chain 
with nearest and next neighbor interactions. N = 2 23 . 

vertex (a;, 1) flips up, given that spins at vertex (x, 2) and level 
x+1 are kept down. Now we relax the spins at levels — 1. The 
spin at vertex (a;— 1, 2) can not flip up in this further relaxation 
unless the spin at vertex (x — 1,1) flips up first. Let Q^ l {h) 
be the probability that the spin at vertex (x — 1,1) flips up, 
but it can not cause the spin at vertex (x — 1, 2) to flip up. We 
denote the probability of both the spins at level x — 1 flipping 
up by Q4 _1 (/i). Therefore, the spins at level x — 1 remain 
down with probability ~ 1 ( - 1 , - 1 ) - ~ 1 (ft) - 1 (ft) . 

We consider the probabilities [P^(l,l), P%(1,— 1), 
Pg(-1, 1) - Qf (ft). -1) - Q%{h), Qf(h), Q x 2 {h) - 



Q%(h) - Q%{h), Q%{h), Q%(h)} as a column vector V x (h), 
so that the recursion relations of the probabilities can be rep- 
resented in matrix form as in Eq. |4] with a 8 x 8 stochastic 
matrix Wh- We find the eigenvector V*(h) for the matrix Wh 
and using V*(h) and p m (h)'s calculate the average magneti- 
zation far away from the boundary. Figure |4] shows the com- 
parison between theoretical and simulation results for Gaus- 
sian quenched random fields with a = 0.5J and a — 1.5J. 
The simulation results are obtained for single realizations, on 
a linear chain of length N = 2 23 with nearest and next neigh- 
bor interactions, with periodic boundary condition. 

In summary, we have demonstrated that there is no jump 
discontinuity in the magnetizations for any quasi-one dimen- 
sional RFIM for unbounded continuous distributions of ran- 
dom fields. We showed that for the RFIM on an n-leg lad- 
der and in ID with interactions extended to n closest neigh- 
bors, for finite n the relaxation of spins at an external field h 
can be described by a stochastic matrix Wh', and the average 
magnetization can be obtained using the eigenvector of Wh 
corresponding to the eigenvalue one, which is a continuous 
and differentiable function of h, for unbounded continuous 
distributions of random fields. We explicitly calculated the 
magnetizations for two simpler cases: (a) the two-leg ladder 
and (b) in ID with nearest and next nearest neighbors inter- 
actions; and confirmed our results using numerical simulation 
for Gaussian distribution of random fields. The question of 
how to take the n — > oo limit, where the magnetizations show 
jump discontinuity below a critical disorder as the external 
field is varied, remains open. 
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